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We present branching fraction measurements of the decays B + — > a+(1260) K° and 
B° -> a^(1260)K + with af (1260) -> tt^tt*. The data sample corresponds to 383 x 10 6 BB 
pairs produced in e + e~ annihilation through the T(4S) resonance. We measure the products of the 
branching fractions B(B + -> a+ (1260) K°) B(a+ (1260) -> 7r-7r+7r+)=(17.4 ± 2.5 ± 2.2) x 10~ 6 
and B(B° ^ a- (1260) K + )B{a- (1260)^ tt + tv--k-)= (8.2 ± 1.5 ± 1.2) x 10~ 6 . We also 



4 



measure the charge asymmetries A c h(B + -> a+(1260) K°) = 0.12 ± 0.11 ± 0.02 and 
Ach{B° -> a^(1260) K + ) = -0.16 ± 0.12 ± 0.01. The first uncertainty quoted is statistical 
and the second is systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Recently the BABAR Collaboration has reported the 
measurement of the branching fraction [l[ and time- 
dependent CP-violation parameters for the process 
5° af (1260) tt t 0]. If this process were mediated by 
a single tree amplitude these measured parameters would 
enable a determination of the angle a of the unitary tri- 
angle of the Cabibbo-Kobayashi-Maskawa (CKM) quark- 
mixing matrix ||. However, in the presence of a penguin 
amplitude with a different weak phase from the tree am- 
plitude the measured angle, called a e g , would differ from 
a [3| . An upper bound on the difference Aa = \a — a c ff | 
can be calculated using SU(3) together with measure- 
ments of the CP-averaged decay rates for either the de- 
cays B+ -> a^~(1260) A or B° -> of (1260) K + and the 
decays B -> Ai(1270) tt and B -> Ai(1400) tt ||. Know- 
ing the value of this difference is important in calculating 
bounds on the angle a. 

There are no experimental measurements of the 
branching fractions of the decays B + — > 0^(1260) K° and 
B a — > (1260) K + to this date. Recent theoretical esti- 
mates of these branching fractions have been calculated 
assuming naive factorization for two different values of 
the mixing angle between the two strange P wave axial 
mesons [6j. The estimated branching fractions lie in the 
range (16 - 52) x 10~ 6 0. Comparison between theo- 
retical predictions and measured quantities is useful to 
test the underlying theoretical hypotheses of factoriza- 
tion and B — > ai(1260) transition form factors. 

We present measurements of the branching fraction for 
the decays B+ -> a+(1260) A and P° -> (1260) K+ 
with af( 1260) — > 7r =F 7r ± 7T ± |||. We also search for a direct 
CP violation by measuring the charge asymmetry A c h, 
defined as (r~ - r+)/(r _ + T+), in the decay rates T± 
for a charged B meson, or T + ( B° — ► a^(1260) K + ) and 
its charge conjugate for a neutral B meson. 

The data were collected with the BABAR detector [9J 
at the PEP-II asymmetric energy e + e~ collider [1QJ. 
An integrated luminosity of 347 fb _1 , corresponding 
to 382.9 ± 4.2 million BB pairs, was recorded at the 
T(45) resonance (on-resonance, center-of-mass energy 
■s/s = 10.58 GeV). An additional 37 fb _1 , recorded about 
40 MeV below the T(4S) resonance (off-resonance), is 
used for continuum background studies. 

Charged particles are detected and their momenta 
measured by a combination of a silicon vertex tracker, 
consisting of five layers of double-sided silicon detec- 
tors, and a 40-layer central drift chamber, both operat- 
ing inside the 1.5-T magnetic field of a superconducting 
solenoid. The tracking system covers 92% of the solid 
angle in the center-of-mass frame. 



Photons and electrons are detected with a CsI(Tl) elec- 
tromagnetic calorimeter. Charged-particle identification 
(PID) is provided by the average energy loss (dE/dx) 
measured in the tracking devices and by an internally 
reflecting ring-imaging Cherenkov detector (DIRC) cov- 
ering the central region. A A"/7r separation of more than 
four standard deviations (a) is achieved for momenta be- 
low 3 GeV/c, decreasing to 2.5 a at the highest momenta 
of the B decay products. 

Monte Carlo (MC) simulations of the signal decay 
modes, continuum, BB backgrounds and detector re- 



sponse [ll| are used to establish the event selection cri- 
teria. Exclusive MC signal events are simulated as B — > 
d (1260) A with a x (1260) -> pw. For the ai(1260) meson 
parameters we take the mass Too = 1230 MeV/c 2 and 
T = 400 MeV/c 2 fillip. We account for the uncertain- 
ties of these resonance parameters in the determination 
of systematic uncertainties. The a^~(1260) — > 7r~7r + 7r + 
decay proceeds mainly through the intermediate states 
(7T7r)p7r and (7r7r)a-7r 13]. No attempt is made to separate 
the contributions of the dominant P-wave (7T7r) p from the 
S-wave (ttw)^ in the channel 7r + 7r~. A systematic uncer- 
tainty is estimated due to the difference in the selection 
efficiency. 

We reconstruct the decay a^~(1260) — * n~ir + n + with 
the following requirement on the invariant mass: 0.87 < 
m ai < 1.5 GeV/c 2 for B° -> of (1260) K+ and 0.87 < 
m ai < 1.8 GeV/c 2 for B+ -> a+(1260) K° . The dif- 
ferent ai mass selections are motivated by charm back- 
ground studies. The intermediate 7r + 7r~ state is re- 
constructed with an invariant mass between 0.51 and 
1.1 GeV/c 2 . Secondary ai(1260) daughter pions are re- 
jected if their PID signatures satisfy requirements for be- 
ing consistent with protons, electrons, or kaons. PID re- 
quirements ensure the identity of the primary charged 
kaon. Candidate K® — > tt + it~ decays are formed from 
pairs of oppositely charged tracks with 0.486 < < 
0.510 GeV/c 2 , having a decay vertex \ 2 probability 
greater than 0.001, and a reconstructed decay length 
larger than three times its uncertainty. 

We reconstruct the P-meson candidate by combin- 
ing an ai(1260) candidate and a charged or neu- 
tral kaon. A P-meson candidate is characterized 
kinematically by the energy-substituted mass toes = 
\J (s/2 + p • p b ) 2 /Pq - p 2 B and energy difference AE = 
Eg — \\fs, where the subscripts and P refer to the 
T(45) and the P candidate in the laboratory frame, re- 
spectively, and the asterisk denotes the T(45) frame. 
The resolutions in toes and AP are about 3.0 MeV/c 2 
and 20 MeV respectively. We require |AP| < 0.1 GeV 



5 



and 5.25 < m ES < 5.29 GeV/c 2 . To reduce fake B- 
meson candidates we require a B vertex x 2 probability 
larger than 0.01. The cosine of the angle between the di- 
rection of the 7r meson from ai(1260) — > pit with respect 
to the flight direction of the B in the ai(1260) meson 
rest frame is required to be between —0.85 and 0.85 to 
suppress combinatorial background. The distribution of 
this variable is flat for signal and peaks near ±1 for this 
background. 

To reject continuum background, we use the angle 9t 
between the thrust axis of the B candidate and that of 
the rest of the tracks and neutral clusters in the event, 
calculated in the center-of-mass frame. The distribu- 
tion of cos Qt is sharply peaked near ± 1 for combinations 
drawn from jet-like qq pairs and is nearly uniform for the 
isotropic -B-meson decays; we require |cos#t| < 0.65. 
The remaining continuum background is modeled from 
off-resonance data. 

We use MC simulations of B°B° and B + B~ de- 
cays to study BB backgrounds, which can come from 
both charmless and charmed decays. The modes 
B° -> 0^(1320) A+and B+ -> a+(1320) A decay to the 
same final states as the signal modes. We suppress these 
backgrounds with the angular variable TC, defined as the 
cosine of the angle between the normal to the plane of 
the 37r resonance and the flight direction of the A me- 
son evaluated in the 3ir resonance rest frame. Since the 
ai(1260) and a2(1320) have spins of 1 and 2 respectively, 
the distributions of the variable Ti for these two reso- 
nances differ. We require \Ti\ < 0.62. 

We have on average 1.3 candidates per event for both 
signal decay modes and we select the B candidate with 
the highest B vertex probability. From the MC simula- 
tion we find that the best candidate selection algorithm 
finds the correct-combination candidate in 92 % of both 
signal decay modes and that it induces negligible bias. 

We use unbinned, multivariate maximum-likelihood 
(ML) fits to measure the yields of B+ -> af(l260)K° 
and B° — ► a|f(1260) A + . The likelihood function incor- 
porates five variables. We describe the B decay kine- 
matics with the two above-mentioned variables AE and 
toeS: as we ll as the invariant mass of the 37r system, a 
Fisher discriminant T and the variable Ti.. The Fisher 
discriminant combines four variables: the angles with re- 
spect to the beam axis in the T(4S I ) frame of the B mo- 
mentum and B thrust axes and the zeroth and second 
angular moments Ln. 2 with respect to the thrust axis of 
the .B-candidate [14J. Since the correlation between the 
observables in the selected data and in MC signal events 
is small, we take the probability density function (PDF) 
for each event to be a product of the PDFs for the in- 
dividual observables. Corrections for the effects of pos- 
sible correlations are made on the basis of MC studies 
described later. The selected data samples besides the 
signal events contain continuum qq and BB combina- 
torial background. 



The BB background has the following components 
in the likelihood: charmless, charm and (1320) A + . 
There are also three additional components: /oA, p°K 
with their yields fixed to the value determined from the 
measured branching fractions [l5[ , and the non-resonant 
p°ir + K with a yield fixed in the fit to the value expected 
using an assumed branching fraction of (2.0±2.0) x 10~ 6 . 
We account for the uncertainties of these branching frac- 
tions in the determination of the systematic uncertain- 
ties. A charged particle from a signal event may be 
exchanged with a charged particle from the rest of the 
event. These so-called self cross feed (SCF) events are 
considered background events. The charmless BB back- 
ground has a dependence on the ML fit observables that 
is similar to that for SCF events, and thus the SCF events 
can be modeled as part of the charmless component. 

The likelihood function is defined as 

N 

C = exp(-^n fc J nX! n i x 0) 

k i=l j 

V k (m ES *)V k (AE i )V k {F)V k (r< )V k (W ) , 

where N is the total number of events in the fit sample, 
n k is the yield fitted for the likelihood component k and 
P k {x l ) is the PDF for observable x in event i . We de- 
termine the PDFs for signal and BB backgrounds from 
MC distributions in each observable. For the continuum 
background we establish the functional forms and initial 
parameter values of the PDFs with off-resonance data. 
The PDF of the invariant mass of the ai(1260) meson 
in signal events is parameterized as a relativistic Breit- 
Wigner lineshape with a mass-dependent width which 
takes into account the effect of the mass-dependent p 
width [l6| . We fix the oi(1260) meson parameters to the 
values found in the branching fraction measurement of 
B° -> a^n T '^. The PDF of the invariant mass of the 
£12(1320) meson is parameterized by a relativistic Breit- 
Wigner distribution. The mES and AE distributions for 
signal are parameterized as double Gaussian functions. 
The AE distribution for continuum background is pa- 
rameterized by a linear function. The tties distribution 
for the combinatorial background is described by an em- 
pirical function that accounts for threshold effects 17 1. 
We model the Fisher distribution T using a Gaussian 
function with different widths above and below the mean. 
The TL distributions are modeled using polynomials. 

In the fit for the decay B° -> <zj~(1260) K+ 
(B + — > af (1260) A' ) there are respectively four- 
teen (twelve) free parameters: five (five) yields and 
nine (seven) parameters affecting the shape of the 
combinatorial background. Table I lists the results of 
the fits. We measure the signal yield bias by generating 
and fitting MC simulated samples containing signal and 
background populations expected from data. The signal 
reconstruction efficiency is obtained from the fraction 
of correctly reconstructed signal MC events passing the 
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selection criteria. Branching fractions for each decay are 
computed by subtracting the fit bias from the measured 
yield, and dividing the result by the efficiency, the 
daughter branching fraction product, and the number of 
BB pairs produced. Equal production rates to 
and B + B~ pairs are assumed. The significance is taken 
as the square root of the difference between the value of 
— 21n£ (with systematic uncertainties included) for zero 
signal and the value at its minimum. 



TABLE I: Number of events TV in the sample, fitted sig- 
nal yield and measured bias (to be subtracted from the sig- 
nal yield) in events (ev.), detection efficiency (e), daugh- 
ter branching fraction product FJ Bi , significance (S) (sys- 
tematic uncertainties included), the products of the branch- 
ing fractions B(B° -> a" (1260) K+) B(a~ (1 260) -> tt+tt-tt-) 
and B(B+ -> o+(l 260) K°) B(a+ (1260) -> 71-^+71-+) respec- 
tively, and charge asymmetry with statistical and systematic 
error. 
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FIG. 2: Projections of a) AE, b) rriES, c) m ai> and d) T for 
the B° — > (17/(1260) K + decay mode. Points represent data, 
dashed lines the continuum and BB backgrounds, and solid 
lines the full fit function. These plots are made with a cut on 
the signal likelihood which includes 30% — 40% of the signal. 
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a+(1260) K" 


N (ev.) 
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9468 


Signal yield (ev.) 
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UBi (%) 
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S(a) 
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6.2 


B(xHT 6 ) 
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FIG. 1: Projections of a) AE, b) oies, c) m ai , and d) T 
for the B + —> af (1260) K° decay mode. Points represent on- 
resonance data, dashed lines the continuum and BB back- 
grounds, and solid lines the full fit function. These plots 
are made with a cut on the signal likelihood which includes 
30% - 40% of the signal. 



In Figs. [T] and [2] we show the AE, iti-es, Tn ai , and F 
projections made by selecting events with a signal like- 
lihood (computed without the variable shown in the fig- 
ure) exceeding a threshold that optimizes the expected 
sensitivity. 



Most of the systematic errors on the signal yield arising 
from uncertainties in the values of the PDF parameters 
are already incorporated into the overall statistical error, 
since they are floated in the fit. The systematic error on 
the B° -> ar(1260) K+ (B+ -> a+(1260) K°) fit yield is 
28 (23) events, which is obtained by varying the PDF 
parameters within their uncertainties. We estimate the 
uncertainty arising from the number of BB pairs to be 
1.1%. The uncertainty in the fit bias correction is 12 
(9) events, taken as half of the fit bias correction. The 
uncertainty in the ai(1260) meson parameters is 12 (6) 
events. The systematic uncertainty assigned to the fixed 
yields in the fit is 3 (4) events. The systematic effect due 
to differences between data and MC for the cos 6>t selec- 
tion is 1.8%. A systematic uncertainty of 2.0 (2.5) % is 
evaluated for the PID. The tracking efficiency contributes 
to the systematics with 1.8 (1.3) %. A systematic uncer- 
tainty of 2.5 % is estimated for the difference in selection 
efficiency in the decay modes through the dominant P- 
wave (7T7r)p and the S-wave (7r7r) CT . The contribution of 
interference between 02 (1320) and ai(1260) is negligible. 
In fact, varying the a2(1320)i^ background with differ- 
ent selection criteria on the angular variable 7i gives no 
significant change to the efficiency-corrected signal yield 
of ai(1260)if. We find also that the systematic effect 
due to different form factors in MC signal simulation is 
negligible. The total systematic error on the branching 
fraction of the charged (neutral) mode is 14% (13%). 

The primary sources of systematic uncertainties in the 
charge asymmetry measurement are the track reconstruc- 
tion or particle identification, the imperfect modelling of 
the interactions with material in the detector and the 
BB background. We study these systematic uncertain- 
ties with MC signal events, qq background in the data, 
and control samples. We assign a systematic uncertainty 



7 



of 0.02 (0.01) to the charge asymmetry of the charged 
(neutral) mode. 

In summary, we have measured the branching frac- 
tions B(B+ -> a+(1260) K°)B(at (1260) -> 7r-7r+7r+) 
= (17.4 ± 2.5 ± 2.2) x 10" 6 and B(B° -> oj" (1260) 
B(a7/(1260) -> tt+tt-tt-) = (8.2 ±1.5 ±1.2) X 10" 6 . 
The charged (neutral) S decay mode is observed 
with a significance of 6.2 (5.1) standard deviations, 
which includes systematic uncertainties. We find no 
evidence for a direct CP-violating asymmetry in these 
decay modes. Assuming B(af(12Q0) — > 7r =F 7r ± 7r ± ) 
is equal to S(af(1260) -> ttMtt ), and that 
B(af(1260) -> (3tt) ± ) is equal to 100% [H, we 
obtain £(5° -> £^(1260)^+) = (16.3 ± 2.9 ± 2.3) x 10~ 6 
and B{B+ -> a+(1260)A° = (34.9 ± 5.0 ± 4.4) x 10~ 6 . 
These results are in reasonable agreement with the 
theoretical estimates. 
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